Gluon condensation and scaling exponents for the propagators in Yang-Mills theory 
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We investigate the infrared (strong-coupling) regime of SU(A'')-Yang-Mills theory on a self-dual 
background. We present an evaluation of the full effective potential for the field strength invariant 
FfivF'^" from non-perturbative gauge correlation functions and find a non-trivial minimum corre- 
sponding to the existence of a dimension four gluon condensate in the vacuum. We also relate 
the infrared asymptotic form of the /? function of the running background-gauge coupling to the 
asymptotic behavior of Landau-gauge gluon and ghost propagators. Consistency between both 
gauges in the infrared imposes a new upper bound on the infrared exponents of the propagators. 
For the scaling solution, this bound reads Kc < 23/38 which, together with Zwanziger's horizon 
condition Kc > 1/2, defines a rather narrow window for this critical exponent. Current estimates 
from functional methods indeed satisfy these bounds. 



I. INTRODUCTION 

The understanding of low-energy QCD has continu- 
ously been advancing over the recent decades. Since a 
variety of methods, gauges and pictures have been used 
to investigate QCD in detail, progress has never been 
Unear but many scenarios have been developed abreast. 
Particularly for the challenging problem of confinement, 
a number of successful scenarios exist in parallel. Con- 
temporary research therefore has to clarify whether these 
scenarios mutually support, exclude or simply coexist 
with one another. In this work, we pursue this ques- 
tion for two scenarios of pure Yang-Mills gauge theory 
based on two different gauges: the Landau gauge and 
the background gauge. 

In the Landau gauge, information about the strongly 
coupled sector and confinement is drawn from low- 
order correlation functions which have been studied 
with functional methods such as Dyson-Schwinger equa- 
tions (DSE) [H-lHl, the functional Renormalization Group 
(FRG) , and stochastic quantization d, ; as well 
as with lattice gauge theory [l^, [HI. Many results 
support the K ugo-Oji ma, or Gribov-Zwanziger confine- 
ment scenarios |12h14| | where the impossibility of colored 
asymptotic states imposes conditions on the infrared (IR) 
behavior of ghost and gluon propagators (see section HII)) . 
Reviews on the different approaches and numerical re- 
sults can be found in [l5|, [l^; for further work see 

The background gauge is particularly useful for com- 
puting the full quantum effective action in an expansion 
in terms of gauge-invariant operators. As the quantum 
effective action governs the dynamics of field expecta- 
tion values, information about confinement is contained 
in the resulting nonlinear and nonlocal quantum equa- 
tions of motion. Simple scenarios for the full effective 
action such as the leading-log model [l9M21^ or more so- 
phisticated dielectric confinement models [22 24] indeed 
find confining potentials for the solution of the field equa- 
tions for static color sources. 



Even though the Landau gauge and the background 
gauge have formally much in common, results from one 
gauge generally cannot easily be transferred to the other. 
The background gauge reduces to covariant Lorenz-type 
gauges including the Landau gauge in the limit of zero 
background. However, most results for the background 
gauge are obtained for zero fluctuation field instead, so 
that direct Landau-gauge information is lost. In turn, 
the full background-field dependence in the background 
gauge cannot be reconstructed solely from Landau-gauge 
correlation functions. For special cases, however, such 
a reconstruction in certain parametric limits is indeed 
possible. This has recently been shown for the case of 
the Polyakov loop [2^ , where Landau-gauge propagators 
have been used to determine a variant of the background 
Polyakov-loop potential. These results compare well with 
a related study in Polyakov-gauge, [26j . 

The resulting connection between the two gauges work 
in both ways [25| : given the Landau-gauge propagators, 
an estimate of the deconfinement phase transition tem- 
perature has been computed in excellent agreement with 
lattice simulations. The other way round, a confining 
Polyakov-loop potential in the background gauge imposes 
a new confinement criterion on gluon and ghost propaga- 
tors in the Landau gauge. This criterion is indeed found 
to be satisfied by these propagators. 

In the present work, we further explore the connec- 
tion between the two gauges, concentrating on the run- 
ning coupling and the effective potential for the square 
of the Yang-Mills field-strength. In particular we focus 
on the interrelation between the asymptotic behavior of 
the gluon and ghost propagator for low momenta and the 
/3 function of the background running coupling. Consis- 
tency of the latter with a strongly interacting infrared 
regime already imposes severe restrictions on the scaling 
exponents of the propagators. Furthermore, we report 
on a non-perturbative calculation of the full effective po- 
tential without any polynomial truncation. This allows 
to search for a non-trivial minimum, indicating the con- 
densation of gluons in the vacuum. 
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This work is structured as foUows: In Sect.lTTl we intro- 
duce a non-perturbative Renorniahzation Group equa- 
tion, the Wetterich equation for the case of non-zero 
background and elaborate on the extension of Landau- 
gauge propagators to the background gauge. We then 
specialize to the case of a self-dual background in 
Sect, mil which allows to evaluate the asymptotic form 
of the /? function of the running coupling. This can be 
related to the asymptotic behavior of the propagators, 
see Sect. ITV)) . Finally, we report on our results for a nu- 
merical evaluation of the full effective potential in SecfVl 
Conclusions and discussion are summarized in Sect. IVII 

II. BACKGROUND-FIELD FLOW 

Information about correlation functions and effective 
potentials is encoded in the effective action F, being the 
generating functional of IPI Green's functions. Nonper- 
turbative access to the effective action is given by the 
functional Renormalization Group (RG) in which the full 
F is constructed successively by integrating out quantum 
fluctuations momentum shell by momentum shell. For 
a generic quantum field theory, this procedure leads to 
a flow equation for a scale-dependent effective action F^. 
Here, k denotes a momentum-scale above which all quan- 
tum fluctuations have been integrated out. The depen- 
dence of Ffe on this momentum scale is determined by the 
Wetterich equation , being an exact functional differ- 
ential equation. Reviews on this topic can be found in 
[28j , for gauge theories see [H, [s^l • The application of 
the background-field method to flow equations in gravity 
is reviewed in e.g. (3]| . 

In gauge theories, momentum-shell integrations in con- 
tinuum quantum field theory can most conveniently 
be performed in gauge-fixed formulations. An effec- 
tive action constructed from only gauge-invariant build- 
ing blocks can nevertheless be obtained by use of the 
background-field formalism [s^, HI] . Here the full gauge 
field is split into a background and a fluctua- 
tion field according to Ap, — + a^. The gauge of 
the fluctuations is then fixed with respect to the back- 
ground field, whereas the whole action remains invari- 
ant under an auxiliary gauge transformation of the full 
gauge held Afj, and the background held A^j,. Let us 
stress that this transformation is not the physical, but 
just an auxiliary gauge transformation. It allows, how- 
ever, to retain invariance of the standard effective action 
T[A] = T[a = 0, A] under physical gauge transformations. 
Within the functional RG setting a = is allowed only 
once all fluctuations have been integrated out, hence the 
price to be paid for the construction of a gauge-invariant 
effective action is its dependence on two gauge fields, 

dilli]. 

The Wetterich flow equation in background-field for- 
malism reads [s^ [sj] : 

dtrk[a,A\ = ^STi- (r^^^°'> [a, A] + RkY' dtRk, (1) 



where F("'"'[a,A] = ^^T,[a,A] and 9, ^ k^. 
The action F^ is a functional which interpolates between 
the microscopic action at the UV cutoff A and the 
full quantum effective action F, i.e., Tk^\ — )■ Sa and 
rfc-i-o r. The solution to the flow equation ^ pro- 
vides an RG trajectory of action functionals F^ that in- 
terconnects the microscopic action with the full effective 
action. The quantity Rk is a regulator function depend- 
ing on an infrared (IR) cutoff k that suppresses propa- 
gation of momenta smaller than k. The trace (STr) runs 
over all internal indices, momenta and fleld components, 
i.e., gluon and ghost degrees of freedom, including a neg- 
ative sign for the ghosts. 

As we will eventually be interested in the background- 
fleld action 

Tk[A]=rk[a = 0,A], (2) 

simpliflcations arise from the fact that only the fluctua- 
tion fleld propagators in a background, being the inverse 
of F^^'"'' [0, ^] , are required on the right-hand side in order 
to extract the flow of Tk[a = 0, A] on the left-hand side. 
In the following we will concentrate on determining these 
inverse propagators in the background gauge with gauge 
parameter a = 0, i.e., in the so-called Landau-DeWitt 
gauge. These are directly related to a covariantization 
of the well-studied (inverse) Landau-gauge propagators 

-'^ Landau (^'^)' which is a well-suitcd gauge for functional 
RG calculations as it is a fixed point of the RG flow 
[3^. Moreover, FRG flows for correlation functions in 
the Landau-DeWitt gauge directly relate to those ob- 
tained within fully g aug e- invariant flows for the geomet- 
rical effective action [33,li3l- It is interesting to note that 
in such a setting the difference between background and 
fluctuation correlation functions is controlled by Nielsen 
identities [1^ , and indeed supports the above relation be- 
tween Landau-gauge propagators and background gauge 
propagators. 

The relation between the inverse propagators in these 
two gauges can be parameterized by 

Ff °)[0,A] =Fil„,,j2?[A]) (3) 

where VIA] is a suitable background covariant differential 
operator, reducing to the Laplacian for vanishing back- 
ground fleld. The covariant derivative with respect to the 
background fleld is given by D'/f = d^S"^ + gf°^^''A1. The 
function /^^ occurring in combination with the Yang- 
Mills field strength in the second term cannot be de- 
termined from the knowledge of the Landau-gauge prop- 
agators alone. It is nonsingular for vanishing argument 
in order to guarantee the proper Landau-gauge limit. For 
backgrounds with vanishing field strength (such as a pure 
Polyakov loop as studied in 25]), the background prop- 
agators are exactly determined. 

In this work, we will particularly be interested in back- 
grounds of nonzero field strengths. For this, we ap- 
proximate the / terms by a minimal reconstruction in 
terms of replacing the dependence on the momentum 
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of the Landau-gauge propagators by the correspond- 
ing background-covariant Laplacians for the transversal, 
longitudinal and ghost modes, -> I?t, I'l, I'gh, respec- 
tively, 

VT^.^ = -D'^Sf,^ + 2z(?F^^, = -D^D^, (4) 

and Pgh = —D^. Here the transversal Laplacian also 
contains the spin-1 coupling to the background field. 
This leads to the following construction of the back- 
ground field inverse propagators 

rfLndaub') -> rr^°)[o, A] = ril„,,j2?), (5) 

where V = I?T,L,gh, respectively. 

For the flow equation ([T]), also the regulator needs to 
be specified. Background-gauge invariance requires a de- 
pendence of the regulator on the background-covariant 
momentum. RG invariance for non-vanishing IR cut-off, 
[29I [35I I , and spectral considerations [s^ , then lead to 

i?, = rf^°)(fc^)r(y), 2/ = ^. (6) 

The regulator shape function r{y) encodes the detailed 
prescription of the momentum-shell integration and will 
be chosen as r{y) = here. 

It has been shown that the form of the regulator can 
be used to maximize the physics content within a given 
approximation pol l4l| . very similarly to the construc- 
tion of improved and perfect actions on the lattice. For 
the numerical computations in the present paper we use 
results for Landau gauge propagators li^ , that have 
been obtained within functional optimization (29j . For 
the general considerations we leave the shape function 
r(y) in ([5]) unspecified. 

Finally, it should be emphasized that the use of back- 
ground field flows is not limited to the pure Yang-Mills 
sector, but has also successfully been applied in QCD 
calculations including quark fluctuations p3l - l46j . 

III. EFFECTIVE-ACTION FLOW IN A 
SELFDUAL BACKGROUND 

We are interested in computing the flow of F^; in a 
derivative expansion, keeping all operators of the type 
{F^yF^u)^ , with n G N. This can be summarized in 
terms of an effective potential W{F'^) which is an an- 
alytic function of F^ — F^lvF^iy For this computation, 
it suffices to evaluate the fiow equation in a background 
of covariantly constant field strength that allows for a 
unique identification of (i^^^-F^i^)" operators. 

To understand why the flow is independent of the spe- 
cific choice of background, consider theory space, which 
is the space of all couplings of operators of a given field 
content, compatible with the chosen symmetries. The 
flow equation defines a vector field in this space by de- 
termining a /? function for each running coupling. After 



a basis is chosen in this space, i.e. the operators are spec- 
ified uniquely as functions of the field content, a specific 
field configuration will only serve to distinguish different 
operators. The (3 function of a specific running coupling 
can be evaluated by projecting the flow equation onto 
a suitable field configuration, but it will not depend on 
this configuration as long as it uniquely determines the 
desired operator. 

In principle, a covariantly constant colormagnetic field 
would be sufficient to project onto the effective poten- 
tial of the field strength invariant. However, the back- 
ground has a second meaning in the present formalism: 
the flow is regularized with respect to the Laplace-type 
spectra in the given background, see Eq. Purely 
magnetic backgrounds now suffer from the problem of the 
tachyonic Nielsen-Olesen mode in the spectrum of fluctu- 
ations [13] typically spoiling perturbative computations 
[48| . The tachyonic mode indicates the instability of the 
colormagnetic background, which may be quantified by 
an imaginary part of the effective action as a decay rate. 
The functional RG can indeed deal with this tachyonic 
mode problem owing to its well-controlled IR regulariza- 
tion [33, Still, the technical complications are sub- 
stantial and the physics of the tachyonic mode in terms 
of, e.g., decay rates of constant magnetic backgrounds is 
not of primary interest. 

Instead, we propose to use the only known stable con- 
stant background field, a covariantly-constant selfdual 
background with F^^ = F^^,, as first analyzed in [i^ 
and [sol [Slj . The covariant Laplacian I?t has nontriv- 
ial zero modes ao ("chromons"), as discussed in detail 
in [50l . [sij . which carry important perturbative contri- 
butions, such as a significant contribution to the 1-loop 
/3 function. As a disadvantage, the selfdual background 
will not uniquely project onto the desired operators, since 
and F'^^,F^jj become indistinguishable. Whereas 
this can potentially lead to systematic errors in the deter- 
mination of the effective potential W(F^), our important 
conclusions drawn from the /3 function of the running 
coupling will be unaffected, see below. 

The covariantly constant selfdual background is given 
by the gauge potential 

A''^ = -^^F^j^^Xyu"- , n°- ^ const., = 1, (7) 

where 

F12 =F^i = f = const. (8) 

Apart from antisymmetric partner components, all other 
components are zero and the field strength is of abelian 
type as the commutator of two gauge potentials vanishes. 
Here, n° is a constant vector in color space that can be 
rotated into the Cartan subalgebra. For later use, we 
define vi (with I = l,...,iV^ — 1) as the eigenvalues of 
^afrpa-^bc _ I^jIJ, Por the traces, we need the 

spectral properties of the relevant operators 2?t, —D^ in 
4 dimensions, which can be found in appendix 1X1 
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In the following, we treat the zero mode ao^, satisfying 
(l?T)^^^aQ^ = 0, separately, as it requires special atten- 
tion for the construction of a suitable regulator. The 
gauge field is decomposed according to 

= A,, + GTfj. + flL,, = + a^^ + aL^ + aofj,, (9) 

where is the selfdual background field, axx are 
transversal and longitudinal gluon modes, and aJp de- 
notes the transversal modes except for the zero mode. 
Our truncation can then be written as 



being the value for the optimized regulator Q- The reg- 
ulator dependence in functional RG cornputations leads 
to a Kc range of Hc € [0.539 , 0.595], see 0; for a specific 
flow, see |8|. In Fig. [TJ we show the momentum depen- 
dence of the ghost- and gluon propagator as obtained 
from a functional RG study jlj in comparison to lattice 
results 

The scaling solution satisfies both, the Kugo-Ojima 
confinement criterion Kc > featuring an infrared en- 
hanced ghost [l^ . as well as the Gribov-Zwanziger sce- 
nario [ij, [13] • The latter relates color confinement to 
horizon conditions for gluons and ghosts which for the 
scaling solution is satisfied if 
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(16) 



(10) 



where c, c denote the ghost fields. 

The inverse propagators for gluons and ghosts are 
taken from Landau-gauge calculations. In the deep in- 
frared, a family of solutions exists that is parameterized 
by the infrared boundary condition for the ghost propa- 
gator 0, im . They can be written as 



ri'ii^[o,o](/)=p2z^(/)nT(p)ii, 



rriL[o,o](p^) 



(11) 



The transversal projector satisfies IIt (p) = Sf^^ — ^^Pii ^ 
and the identity is understood to apply to color indices. 
Herein the wave function rcnormalizations 



Za{p' ^ 0) ~ ip') 



2\KA 



Zc{p^ ^ 0) ip'r^ (12) 



are parameterized by critical exponents ka and Kc in the 
deep infrared. In the Landau gauge, the longitudinal 
propagator remains trivial as the longitudinal modes de- 
couple completely from the flow. 

The so-called decoupling solution, found in lattice stud- 
ies of Landau-gauge Yang-Mills theory [ll|, and also in 
continuum studies, (5l.[l7|. is characterized by the critical 
exponents 



-1 



and 



0. 



(13) 



which corresponds to a positivity-violating gluon propa- 
gator [H, [i^l, indicating the confinement of gluons. 

The scaling solution, first identified in [l|, has only 
one independent critical exponent the non- 

renormalization theorem for the ghost-giuon vertex [ssf 
(proven to all-orders in perturbation theory) implies the 
sum rule 

KA = -2kc (14) 

in d = 4 dimensions 0, 13 • most functional computa- 
tions we are led to {d = 4) 



0.59535. 



and KA 



-1.1907..., (15) 



The gauge correlation functions can also be used to re- 
late color confinement to quark confinement: as shown in 
[2^ , the infrared behavior of the propagators determines 
the effective potential of a Polyakov-loop order param- 
eter for the deconfinement phase transition. From the 
behavior of the effective potential, a sufficient criterion 
for quark confinement can be deduced, implying 



1 

Ka>- 



(17) 



for the scaling solution. The decoupling solution also sat- 
isfies a corresponding confinement criterion. Hence both 
the decoupling and the scaling solution for the propaga- 
tors correspond to scenarios where all color charges are 
confined. The scaling solution, however, is the only so- 
lution which is comaptible with global BRST invariance, 
see 0, [sJl ■ For further work on the infrared behaviour of 
Landau-gauge Yang-Mills theory see, e.g., [isj . 

The asymptotic forms of the propagators can be pa- 
rameterized as 



rS^P^) = r^^(p^)HT(;,)ii, r)::>{p^) = n%{p^)n, 

(18) 

where the scalar functions ^''^\/cip'^) are the IR- 
regularized generalizations of Eq. (|TT|) (42j , 



r^^\ip')=lA 



{p^+CAkY^^' 

(A|cd)^^ 



7c 



(19) 



Here, the 7a,c's are simple proportionality constants 
which account for the difference between Aqcd and the 
scale where this asymptotic form takes over. The con- 
stants CA,c arc manifestations of the regulator depen- 
dence in the asymptotic regime; generically, they are of 
order 0(1). In the absence of any IR regularization, i.e., 
fc — > or CA,c 0, Eq. p9|) reduces to the standard 
form, cf. Eqs.' pT|) . P^ . 

For the regularization of the zero modes, some care 
is required, since the seeming standard choice T> = 2?t 
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would imply that the regulator function r{y) is always 
evaluated at the singular point y = V/k^ (because 
I?T = on the zero-mode subspace). In particular, this 
would not properly account for the decoupling of the zero 
mode, once it acquires a field-dependent mass, which al- 
ready happens pcrturbativcly [SOj. Instead, we choose 

I 



V = —D^ as the argument of the regulator on the zero- 
mode subspace which makes the regulator satisfy all stan- 
dard requirements. On the zero-mode subspace, we have 

V = -D'^ 2fi, cf. first line of Eq. for n = m = 0. 
The flow equation for the background potential in the 

selfdual background can finally be written as 



+ hTdtRk,Li-D^) (ri'i + RkA-D'')y^ + ^tr PodtRk.AhD^) {^kii + Rk,Ai-D^)) ' , (20) 



where Pq denotes the projector onto the zero-mode sub- 
space. The zero-mode trace simply corresponds to the 
terms to = n = in the spectrum (jAip . The subscripts 
A, c and L in the regulator function imply the use of the 
appropriate transversal, ghost or longitudinal two-point 
function in its argument. 

The trace measure factor in Eq. ([20]) corresponds to 
the density of states in a selfdual background field, given 
by [.///(27r)]2 [13. 



IV. RELATION BETWEEN THE (3 FUNCTION 
AND THE CRITICAL EXPONENTS 

In the following we will investigate the asymptotic form 
of the /3 function of the background running coupling, 
defined as the prefactor of the F^^^F^^ term: 



(21) 

where n denotes the spacetime volume. For asymptotic 
freedom, the Gaufiian fixed point of the /3 function has 
to be UV attractive (i.e. IR repulsive). On the other 



hand, a minimum requirement for the confinement of 
color charges is an interacting theory in the long-range 
limit. This requirement is incompatible with an IR at- 
tractive Gaufiian fixed point. As we show below, this 
seemingly elementary condition imposes a new confine- 
ment criterion on the exponents ka and Kc- 

According to Eq. (pi]) the (3 function can be evaluated 
from the first term in a Taylor expansion of the effective 
action in powers of the background field 



dtg' 



-g 



1 



'dt-^-^d,r,[A] 



p 



(22) 



Here |y2 denotes the projection onto the Taylor coefficient 
at order Although there can be ambiguities in the 
projection onto (i^^i^F^i/)" and {F^,^F^A'^ for a selfdual 
background, this does not affect the (3 function, as the 
second operator is only non-zero for even powers n due 
to parity conservation. 

For the evaluation of Eq. (P^ . some care is needed as 
the trace over the Laplace-type spectra and the projec- 
tion onto the p order do not commute. As the degener- 

acy factor in the trace already carries a factor of fj^^y , 
all / dependence outside the operator trace can already 
be ignored due to the projection, yielding 
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.(2) 
kA 



(fc2 



n(2) 



(fc2 



9*rg(fc2) / r^^^(o) 



(2), 



(2), 



1 tre" 



.(2) 



(0) 



(2) 



tr e 



(23) 



To obtain the /3 function, we extract the coefficient of 
the expansion of the heat-kernel trace over coordinate 
and color space at second order in /, 



Tr,.r,e 



(4-)^ h Binh^(it) 



(47r)2 ^ V 3 



(2) 

Evaluating dtTj^ with the infrared asymptotic form of 
the inverse propagator Eq. p9p yields the infrared form 
of the 13 function, whereas using the perturbative form 
of the inverse propagators r(2)(p2j ^ ^2 j-eguH^g the 
standard one-loop form of the /3 function. Recovering 
this universal term within our setting can be viewed as a 
simple check of our formalism. 

Focusing on the infrared asymptotic form of the prop- 
agators Eq. (Uni), we obtain to lowest order in the 
background-field coupling 

(47r)2 \^Ck + 1 3 3 Ck + 1 J 

(25) 

where we have used that X^T" ^ '^f ^ -^c- Here we have 
defined 



.(2) 
kA 



(0) 



.(2) 
kA 



(fc2) 



CA 



l + CA 



1 + KJ 



(26) 



Incidentally, the perturbative 1-loop /S function /3g2 = 
~T (ItO^ obtained from Eq. (^5)) in the limit of per- 
turbative critical exponents ka,c — ^ and in the absence 
of an additional IR regularization ca,k 0. 

For retaining an interacting theory in the long-range 
limit, the /? function near the Gauf5ian fixed point must 
not change sign. This gives us a new bound on the IR 
values of the k's, 



-c)-J 



,}_+_KA 

' C.- + 1 



> 0, (27) 



For the accidental case of an equality, the sign of the (3 
function would have to be determined from higher order 
terms in and a corresponding additional bound would 
follow from the prefactor of these higher-order terms. 



Let us first concentrate on the scaling solution: Us- 
ing the sum rule ka = — 2kc, this bound yields a crit- 
ical value of Kc,crit(cA) as a function of the regulator- 
dependent constant ca- For an interacting IR theory, 
the true critical exponent has to satisfy 



t{cA) 



(28) 



such that our criterion imposes for the first time a rele- 
vant upper bound on the Landau-gauge exponents. The 
strongest bound (minimum upper bound) is achieved in 
the limit of ca — > oo, where 



rit(cA OO) 



23 
38 



0.6053. 



(29) 



The maximum upper limit for Kccrit is reached for ca ~ 
0.1073, where Kc.crit — 0.72767. For completeness let 
us mention that, for values of ca ^ 0.2, the inequal- 
ity Eq. (P?]) can also be satisfied if Kc > Kc.crit2, as the 
nonlinear inequality Eq. ((27|) bifurcates. We find that al- 
lowed values for Kc lie to the left of the red/upper curve 
(see right panel of Fig. [2]). For < ca < 0.1073, only 
the bound Kc < 1 from unitarity [H remains. How- 
ever, let us stress that the limit ca corresponds to 
a highly asymmetric regularization as the contribution 
from transverse gluons to the 13 function is removed in 
this limit (note that c„ — > cio for ca — ^ and ka < ^1)- 
For this case, the Gaufiian fixed point is naturally IR 
repulsive for all values of Kc > ^. 

To summarize, our argument based on the mere ex- 
istence of an interacting IR regime imposes a bound on 
the IR critical exponent Kc < Kc.crit- Both sides of this 
inequality are regulator dependent, such that the bound 
has to be satisfied in every regularization scheme as a nec- 
essary criterion. A sufficient criterion is therefore given 
by imposing that Kc < Kc,crit,min = 23/38 ~ 0.6053 holds 
for all regularization schemes. 

This bound is furthermore required to be satisfied 
due to the following argument: an IR fixed point for 
the background-field running coupling is observed in the 
background-field gauge in a derivative expansion f37| as 
well as for the running coupling fixed at the ghost-gluon 
vertex in the Landau gauge in a vertex expansion 
Both running couplings are defined with the aid of non- 
renormalization theorems arising from gauge invariance. 
From the conjecture that both running couplings are ac- 
tually linked on all scales, consistency requires the ex- 
istence of the fixed point in both gauges. Even though 
such a fixed point is not observed in the simplified trun- 
cation leading to Eq. (P^ . it is expected to arise in a 



7 



0.70 
0.68 
0.66 
0.64 
0.62 
0.60 



I.Oh 




0.05 0.10 0.15 0.20 



FIG. 2. ftccrit as a function of ca- Left panel: continu- 
ous branch of Kc.crit that exists for all ca (all regularization 
schemes); for the ca values displayed here, ca ■> 0(1), Kc,crit 
represents an upper bound for the ghost critical exponent Kc- 
Right panel: for 0.1 < ca ^; 0.2, the bound bifurcates, such 
that all Kc values to the left of the red/upper curve are al- 
lowed. The horizontal line marks the value Kc — 0.595 as 
found in must functional calculations for the scaling solution 
which obviously satisfies our confinement criterion. 



larger truncation as used in [37| . But already without 
knowing the full IR flow of the background field coupling, 
consistency of the leading-order result of Eq. (^5]) with 
the Landau-gauge coupling gives rise to the bound on 
the scaling exponent. This consistency would be spoiled 
if the above first coefficient of the /3 function changed 
sign. Therefore, we infer that Kc < || is also required 
for the consistency between the background-gauge and 
the Landau-gauge results, independently of higher-order 
terms that generate an IR fixed point of the coupling. 

Physically our upper bound for may be interpreted 
as a criterion for an IR interacting theory and thus ul- 
timately for confinement. Therefore, we conclude that 
for He larger than some critical value Yang-Mills theory 
would not be confining. Therefore it is very reassuring 
that all results for Kc obtained by functional methods sat- 
isfy the sufficient criterion bound Kc < || , c.f. Eq. ([T5|. 

It is important to note that also that the decoupling 
solution for the propagators = ^ 1, = that has 
been found in many lattice simulations (and can possibly 
be understood as a different way to deal with the gauge 
ambiguity in the Landau gauge, see [s^ ) fulfills our con- 
finement criterion which in this case is given by the more 
general bound Eq. (P7| . 



V. GLUON CONDENSATION FROM THE 
EFFECTIVE POTENTIAL 

Dimensional transmutation and the scale anomaly go- 
ing along with the quantization of Yang-Mills theory 
require a non-zero expectation value for the energy- 
momentum tensor. This implies a non-trivial vacuum 
structure: 



Indeed phenomenological estimates indicate a non-zero 
value [Sa], being interpreted as a condensation of gluons 
in the vacuum. These findings have been corroborated 
by^other methods [13, Ull , as well as lattice gauge theory 

The effective action for a colormagnetic background 
field has been evaluated at one- loop order in [60j . where 
a non-trivial minimum has been found (Savvidy vacuum). 
Due to the tachyonic mode of the gluon propagator on 
such a background this configuration is unstable j47| . 
which can be amended by introducing a spatial inhomo- 
geneity into the field configuration. This vacuum configu- 
ration is referred to as the Copenhagen vacuum [13, l6ll - 
l63j . A more severe problem is posed by the fact that 
the minimum lies in the non-perturbativc domain, i.e. 
beyond the perturbativc Landau pole of the Yang-Mills 
coupling, questioning the relevance of a perturbative es- 
timate. 

First non-perturbative studies from functional meth- 
ods indeed provided indications for gluon condensation 

M- 

Here we will further pursue this question, based on 
the knowledge of full correlation functions. As a simple 
parametrization of these correlation functions, we still 
use the asymptotic form displayed in Eq. (|19p which we 
amend with the fc-dependent critical exponents K,A{k) 
and Kc(fc) in accordance with the propagators in [4^ . A 
suitable interpolation between KA,c{k — > oo) — > and the 
corresponding IR values KA.c{k —J- 0) — > ka.c of Eq. 
can parameterize the full momentum dependence of the 
correlation functions. This allows to evaluate the effec- 
tive potential from an integrated form of the flow equa- 
tion: 



(2) 



(31) 



where fcuv 

p2 _ pa pa 



IS an 

f2 



initial ultraviolet scale, and 
fiu^ — 4/^. The initial condition Wk^^^F"^) is 
fixed deep inside the perturbative regime at fcuv = 10 
GeV, 



uv 



W, 



F' 



fc=10GoV 



167r0.2294' 



(32) 



{F') ^ 0. 



(30) 



where we have used the peak of the maximum of the 
gluon dressing function 1/Za{p^), see left panel of Fig.[Tl 
for relating our YM scales to that used in lattice compu- 
tations. The normalization is such that the related string 
tension a (computed on the lattice) is given by ^/a ~ 440 
MeV. Our initial condition is thus self-consistently de- 
termined from our main input, the Landau-gauge ghost 
and gluon propagators, which goes along with a coupling 
as{k = lOGeV) ~ 0.2294 at the initial scale. At this 
initial cut-off scale, k = lOGeV, higher-order operators 
do not contribute significantly to the effective potential, 
which is determined exclusively by the functional form 
of the classical action. Alternatively, an RG-improvcd 
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initial condition in the form of the one-loop effective po- 
tential could be used. 

Eq. (|3T|) can then directly be integrated numerically. 
For SU(3), we find a non-trivial minimum at 



F'^ w 0.93 GcV'' (3.46AQCD)^ 



(33) 



where we have used Aqcd = 284 MeV, cf. solid line in 
Fig. El 

We emphasize that the formation of the non-trivial 
minimum is mainly driven by the dynamics in the mid- 
momentum regime, the form of the propagators in the 
deep infrared is not crucial here. The existence of a 
gluon condensate in the vacuum thereby occurs for both 
the scaling and the decoupling solution, and its value 
does not significantly depend on the infrared asymp- 
totics. This agrees with the observation in [2^, where 
quark confinement can also be inferred from the propa- 
gators, independent of their infrared asymptotics. 

Note that on the right-hand side of the flow equation 
the color eigenvalues enter. We may choose these 
along one of the two directions of the Cartan subalgebra, 
which yields an uncertainty of about 10% in the value of 
the minimum. 

As a simple parametrization, we fit the effective poten- 
tial to a function of the form 



W{F^) = aF^ ln(&F2), 



(34) 



which is inspired by the corresponding one-loop results 
[60| with two fit parameters a and 6, cf. dashed line in 
Fig. El. The fit yields a = 0.00528 and b = 0.433GeV~^ 



0.010 



0.005 




FIG. 3. The effective potential for SU(3) as a function of 
(thick bhie hne), and the one- loop inspired fit to the numerical 
data of the form aF^ In bF^ (orange dashed line). 

With some reservations, the one-loop result might thus 
be interpreted as providing indeed a reasonable qualita- 
tive estimate of the true nonpertubativc functional form 
of the effective potential. Let us stress once more that 
the one-loop calculation is, of course, not reliable, as the 
predicted minimum lies in a regime inaccessible to per- 
turbation theory. It is only its qualitative prediction of 
the functional form of the potential which seems to be 
altered little by non-perturbative corrections. 



A few comments on the numerical values for the gluon 
condensate are in order: first of all, as the selfdual back- 
ground does not distinguish between the invariants 
and FF, our result for the condensate receives con- 
tributions from both operators; under the assumption 
that condensates of both types exist, the value given 
in Eq. p3p should be considered as an upper bound 
for the phenomenologically more relevant F^ conden- 
sate. In this sense, our result {F'^)/Att ~ 0.074GeV'* is 
well compatible with recent phenomenological estimates, 
(i^2)/47r ~ 0.068(13)GeV'' from spectral sum rules 
(note that our field definition differs from that of [6^1 by 
a rescaling with the coupling, cf. Eq. ((2T|) ). The good 
agreement might even indicate that the contribution from 
FF operators is rather small. Also, we expect the inclu- 
sion of dynamical quark degrees of freedom to decrease 
our condensate value slightly, owing to their screening 
property. 

We have used a background that can only be consid- 
ered as an approximation of the true ground state lo- 
cally. For studies involving a more realistic gauge-field 
configuration than the one considered here see [6J] and 
references therein. 

Within the leading-log model which is based on the as- 
sumption that the infrared effective action is essentially 
given by its one-loop functional form (j34p . the conden- 
sate value is related to the string tension in a simple man- 
ner [l^: the static potential between two opposite color 
charges arising from the nonlinear field equations follow- 
ing from Eq. (|34p grows linearly with distance, V = err, 
where the string tension a and the minimum of the ac- 
tion obey, y/a ~ (Qk)^/*. Here, Q = 4/3 is a simple 
color factor for SU(3), and k = (l/4)F^|,„in denotes the 
minimum of the action, i.e., the condensate value. From 
our result p3p . we obtain the estimate 



.1/2 



747MeV 



(35) 



which, in view of the rather restrictive truncation and 
the simplicity of the leading-log model, compares rather 
favorably with the value cr^/^ ~ 440MeV which we used 
for fixing the scale of the propagators in the first place. 
We should stress that the string-tension a appears in two 
very different meanings in our calculation: It first occurs 
as an input scale for fixing the initial condition for the 
fiow equation, i.e., it fixes the absolute scale of the prop- 
agators. Wc then derive the physical string tension in a 
nontrivial way from the minimum of the effective poten- 
tial via the leading-log model. This output is linked to 
a mechanism of confinement, and has therefore acquired 
a physical meaning beyond pure scale fixing. With re- 
gard to the approximations involved in our calculation 
and, in particular, in the confinement model used to map 
our results onto the physical string tension, the order- 
of-magnitude equivalence of "input scale" and "output 
physics" is satisfactory. 
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VI. CONCLUSIONS 

Based on the knowledge of low-order gauge correla- 
tion functions from Landau-gauge calculations, we have 
used an approximate mapping onto propagators in the 
Landau-DcWitt gauge in a background field. This al- 
lowed us to extract nonperturbative information about 
the effective action of Yang-Mills theories. We have spe- 
cialized to a selfdual constant background which - apart 
from technical simplicity - is known to be stable against 
fluctuations and thus a candidate for at least local ap- 
proximations of the Yang-Mills vacuum. 

We have concentrated on two properties of the effective 
action: the running of the background field coupling and 
the form of the effective potential for the self-dual field 
strength. From the simple requirement that the back- 
ground field coupling should remain finite, and the the- 
ory thus interacting, in the long-range limit, we have de- 
rived a nontrivial criterion for the infrared behavior of 
Landau-gauge correlation functions. As interactions arc 
mandatory for confinement, our criterion may be viewed 
as a confinement criterion. For the scaling solution of 
these correlation functions which is characterized by a 
single critical exponent Kc, this criterion translates into 
a new upper bound on Kc- In its sufficient version, the 
confinement criterion is satisfied if 

23 

< tccrit,min = - 0.6053, (36) 
oo 

where Kc^crit.min is a minimum value obtained in a specific 
regularization scheme. If Kc satisfies this bound in any 
scheme, the system remains interacting in the infrared. 
In general, both sides of the inequality are regularization- 
scheme dependent, such that the condition that Kc[Rk] < 
Kc,crit[-Rfc] is a necessary criterion for confinement which 
then has to be satisfied for each regulator Rk. 

Read together with the confinement criterion for quark 
confinement as derived from the Polyakov loop potential 
in [25j . Kc > J, and the horizon condition for color con- 
finement [3l , Kc > ^ , our new criterion defines a window 
for the infrared critical exponent, 0.5 < Kc < 0.6053, 
thereby tightly fixing the possible infrared behavior of 
the gauge correlation functions for the scaling solution in 
the infrared. It is reassuring that the result Kc ~ 0.595 
obtained from many functional calculations lies precisely 
in this window. 

Furthermore, we numerically investigate the full 
effective potential for the field strength invariant 
p2 _ F^^F^j^ and find a non-trivial minimum at 
F2~0.93GeV^. Our calculation is fully non- 
perturbative, and therefore supports the conclusions 
drawn from the one-loop effective action (60j in a non- 
trivial way. 

Most importantly, effective models of confinement such 
as the leading-log model or dielectric confinement models 
receive strong support from our results. These models 
map the nontrivial vacuum structure onto classical field 



equations. For instance, it can be shown rather straight- 
forwardly in the leading-log model that a nontrivial mini- 
mum of the effective action at nonzero F'^ is already suffi- 
cient to produce a linearly rising potential for static color 
charges at long distances [l^. Our calculation therefore 
provides for an explicit example how different pictures 
of confinement can not only exist in parallel but actu- 
ally support each other, potentially being two sides of 
the same coin. A similar observation has been made for 
the Coulomb gauge in the Hamiltonian approach and the 
Gribov-Zwanziger confinement scenario |67| . 

Quantitatively, our estimate of the selfdual gluon con- 
densate ~ 0.93 GeV^ compares favorably with phe- 
nomenological estimates [65| (F^ ~ 0.85 GcV^). It is 
also in agreement with the large-distance limit of the 
static potential between a quark and antiquark. yielding 
a string tension of cr^/^ ~ 747MeV within the leading-log 
model of confinement [l^ . Our nonperturbative calcula- 
tions, starting from the microscopic theory, thus can be 
interpreted as providing a fundamental justification of 
the two basic assumptions of the leading-log model: the 
functional form of the effective action and the free 
parameter given in terms of the minimum of the effective 
action, i.e., the gluon condensate. Our studies there- 
fore provide an explicit example of how seemingly differ- 
ent confinement scenarios, the Kugo-Ojima and Gribov- 
Zwanziger scenario on the one hand, and the leading-log 
model on the other hand, can not only exist in parallel 
but actually support each other. 

Our results arc affected by several sources of uncertain- 
ties: first of all, the mapping of Landau-gauge propaga- 
tors onto background- field propagators is not unique, but 
has been realized in a minimal-coupling fashion. Whereas 
we have concentrated onto lowest-order correlation func- 
tions in the Landau gauge, also higher orders contribute 
to the full effective potential in the background gauge. 
In principle, it is straightforward to include such higher- 
order correlation functions, as, e.g., computed in [66| . 
into our formalism. Also, differences between gluonic 
field invariants, such as (F^)" and (FF)", which have 
been ignored for the selfdual background can be treated 
with suitable heat-kcrncl expansion techniques in future 
studies. 
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Appendix A: Spectral properties of Laplace- type needed in this work are: 

operators , 

Spccj -D^[ =2fi{n + m+l), n, m = 0, 1, 2, . . . 



Let us summarize the spectral properties of the 
Laplace-type operators involved. We concentrate here on 
the particularities of the selfdual background. The gen- 
eral trace technology in covariantly constant backgrounds 
as it is relevant for the computations in this paper can 
be found in [sj, [s^- The selfdual-background spectra 



c! /t:, \ _ / 2./K'^ + + 2) , multiplicity 2 
spec|i^Tj> y2fi{n + m) , multiplicity 2' 

(Al) 

where /; ~ \i^i\f, and i^i denotes the eigenvalues of the ad- 
joint color matrix n°'T'^. The covariant spin-1 Laplacian 
I?T has a double zero mode for n = to = which is due 
to the symmetry between colorelectric and colormagnetic 
field. Using the trace technology of [13, [13] , the spectral 
problem of the longitudinal Laplacian I'l can be mapped 
onto that of — D^, such that Eq. (|A1[) is sufficient for the 
calculation in the main part of the paper. 

Defining Tr as the trace without the zero mode, we 
make the following useful observation for the trace over 
some function F: 



""^/-cX^foo cjocjo oo 

Tr;,L-^(PT)=2 ^ (f ) E -^(2/;(n + m + 2))+^^ J-(2/,(n + TO))+^^(2/,n) 



1=1 

N^-l 



1=1 



\ n,m— 
2 oo 



n— m— 1 



n=l 



= 4E (^) E HVdn + m + l)) 



n.m— 



(A2) 



where the trace subscripts denote traces over coordinate 
space "x", color space "c" and Lorentz indices "L". In 
other words, there exists an isospectrality relation be- 
tween —D^ and the non-zero eigenvalues of I?t- Similar 
isospectrality properties are also known for the Dirac op- 



erator in a self-dual homogeneous background [68|. As 
a consequence, all gluon and ghost modes except for the 
two zero modes couple in the same fashion to the self- 
dual background. This basic observation has been used 
for the decomposition of terms in Eq. (^0)) . 
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